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ABSTRACT 

We measured the angular diameter of the hthium-rich K giant star HD 148293 using Georgia 
State University's Center for High Angular Resolution Astronomy (CHARA) Array interferome- 
ter. We used our measurement to calculate the star's effective temperature, which allowed us to 
place it on an H-R di agram to compare it with other Li-r ich giants. Its placement supports the 
evidence presented bv ICharbonnel fc Balachandran that it is undergoing a brief stage in 

its evolution where Li is being created. 



Subject headings: infrared: 
individual: HD 148293 

Introduction 



stars, stars: fundamental parameters, techniques: interferometric, stars: 



G and K giant stars are expected to have low 
lithium abundances because Li is destroyed in 
main-sequence stars except in the outermost lay- 
ers (1% - 2% by mass). As the star ascends 
the red giant branch of the H-R diagram, the 
convective envelope deepens, diluting the exist- 
ing Li and furthe r reducing the observable Li 



( Brown et al. |[l98i). This effect is seen in most 



G and K giant stars, though a small minority 
1%) exist that are unexpec te dly rich in Li (e.g., 
Wallerstein fc SnedenI 19821: Pilachowski et al 
jranl ll 



IQQdiFekel fc Balachandranlll993l ). 

Lithium abundance calculations are sensitive to 
temperature variations, so knowing the effective 
temperature (Tcff) of a star is vital in estimat- 



ing its Li co ntent. HD 14 8 293 w as discovered to 
be Li rich bv iBrown et~aLl (|l989l ). who estimated 
its temperature using published photometry and 
color-te mperature relations. Their va l ue wa s also 
used bv lCharbonnel fc BalachandranI ( 2000l here- 
after CBOO), who placed HD 148293 on an H-R 
diagram and found it was at a phase known as the 
"bump in the luminosity function" . This phase is 
characterized by an outwardly-moving hydrogen 
shell, leading to a short-lived phase of Li produc- 
tion before it is rapidly destroyed as evolution con- 
tinues. Only low-mass stars that contain a highly 
degenerate helium core and later experience the 
helium flash pass through this stage and spend a 
mere 3% of their ascent on the red giant branch 
at the bump (-80,000 years, CBOO). 
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By directly measuring the angular diameter of 
HD 148293, we are able to calculate its Tcff when 
combined with other observed quantities, such as 
interstellar absorption and bolometric corrections. 
We then modified the H-R diagram presented in 
CBOO to support their claim of proximity to the 
red-giant bump. Section 2 describes our observing 
procedure, Section 3 discusses how HD 148293's 
angular diameter and Toff were determined, and 
Section 4 explores the physical implications of the 
new measurements. 

2. Interferometric observations 

Interferometric observations were obtained us- 
ing the CHARA Array, a six element Y-shaped 
optical- infrared int erferometer located on Moun t 
Wilson, California ( ten Brummelaar et al. 2005[) . 



All observations used the pupil-plane "CHARA 
Classic" beam combiner in the if'-band at 2.14 /im 
while visible wavelengths (470-800 nm) were used 
for tracking and tip/tilt corrections. The ob- 
serving procedure and data re duction process 



emplo yed here are described in iMcAlister et al 
(|2005l ). We observed HD 148293 over two nights 



using two telescope pairs with different baseline 
lengths: 30 July 2010 using the E2-W2 pair with a 
baseline of approximately 156 m and 31 July 2010 
using the W2-S2 pair at approximately 177 m0 

Two calibrators (HD 145454 and HD 147321) 
were selected to be single single stars with ex- 
pected visibility amplitudes >95% so they were 
nearly unresolved on the baselines used, which 
meant uncertainties in the calibrator's diameter 
did not affect the target's diameter calculation as 
much as if the calibrator star had a significant 
angular size on the sky. We interleaved calibra- 
tor and target star observations so that every tar- 
get was flanked by calibrator observations made as 
close in time as possible, which allowed us to con- 
vert instrumental target and calibrator visibilities 
to calibrated visibilities for the target. 

To check for possible unseen close companions 
that would contaminate our observations, we cre- 



^The three arms of the CHARA Array are denoted by their 
cardinal directions: "S" , "E" , and "W" are south, east, and 
west, respectively. Each arm bears two telescopes, num- 
bered "1" for the telescope farthest from the beam com- 
bining laboratory and "2" for the telescope closer to the 
lab. The "baseline" is the distance between the telescopes. 



ated spectral energy distribution (SED) fits based 
on published UBVRIJHK photometric values 
obtained from the literature for each calibrator 
to establish diameter estimates. We combined 
the photometry with Kurucz model atmosphere^ 
based on Toff and log g values to calculate limb- 
darkened angular diameters for the calibrators. 
The stellar models were fit to observed photom- 
etr y after conv e rting magnitudes t o fluxes us- 
i nglcolina et al.l (|l996i UBVRI) and[Cohenelal 



(|2003i JHK). The photometry. Toff and log g val 



ues, and resulting limb-darkened angular diame- 
ters for the calibrators are listed in Table [1] There 
were no hints of excess emission associated with a 
low-mass stellar companion or circumstellar disk 
in the calibrators' SED fits (see Figure [T]). 



HD 145454 




5000 10000 20000 

VWVELENGTH (Angstroms) 




5000 10000 20000 

WAVELENGTH (Angstroms) 



Fig. 1. — SED fits for the calibrator stars HD 
145454 and HD 147321. The diamonds are fiuxes 
derived from UBVRIJHK photometry (left to 
right) and the solid lines are the Kurucz stellar 
models of the stars. See Table [1] for the values 
used to create the fits. 



^Available to download at http:/ /kurucz. cfa.harvard.edu| 
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Table 1 
Calibrator Information. 



PcircLniGtsr 




JJT) 1 /17Q01 

1 1 IJ 14: / OZ i 


Source 




U magnitude 


5.35 


6.22 


Mermilliod aggi) 




B magnitude 


5.42 


6.07 


Mermilliod agQl) 




V magnitude 


5.44 


5.99 


Mermilliod aggi) 




R magnitude 


5.46 


5.99 


Monet et al. l'2003) 




/ magnitude 


5.50 


5.98 


Monet et al. ('20031 




J magnitude 


5.37 


5.79 


Cutri et al. C2003) 




H magnitude 


5.43 


5.82 


Cutri et al. f2003) 




K magnitude 


5.43 


5.77 


Cutri et al. f2003) 




Teff (K) 


9772 




AUende Prieto & Lambert 


(1999) 


log g (cm s^^) 


4.13 




AUende Prieto & Lambert 


C1999) 


Teff (K) 




8600 


Lafrasse et al. (2010) 




log g (cm s~^) 




4.2 


Lafrasse et al. (2010) 




6Ild (mas) 


0.268 ± 0.015 


0.240 ± 0.010 







3. Determination of angular diameter and 

The observed quantity of an interferometer is 
defined as the visibility {V), which is fit to a model 
of a uniformly-illuminated disk (UD) that repre- 
sents the observed face of the star. Diameter fits 
to V were based upon the UD approximation given 
by y = 2Ji{x)/x, where Ji is the first-order Bessel 
function and x = ttB^udA"^, where B is the pro- 
jected baseline at the star's position, ^ud is the 
apparent UD angular diameter of the star, and 
A is the effective wave length of the observation 
(jShao &: Colavital 119921 ). A more realistic model 
of a star's disk involves limb-darkening (LD), and 
rela tionship incorporating the lin ear LD coefficient 



Hx (jHanburv Brown et al.lll974l ) is 



(l-/iA) 



V 



Ji(x) 



MA 



MA 

3 



1/2 



>^3/2(x) 



(3/2 



(1) 



TablelDlists the Modified Julian Date (MJD), pro- 
jected baseline (B) at the time of observation, pro- 
jected baseline position angle (0), calibrated visi- 
bility (V), and error in V (aV) for HD 148293. 

Th e LD coefficient was obtained from lClaret et al 
(|l995[ ) after adopting the Tcff and log g values for 
HD 148293 and the resulting LD angular diameter 



is listed in Table [31 The errors on the UD and LD 
diameters are 4% and the difference between the 
UD and LD diameters is on the order of a few 
percent, and the final angular diameter is little 
affected by the choice of Additionally, the 

combination of the interferometric measurement 
of the st ar's angu lar diameter plus the Hipparcos 
parallax (jvan Leeuwen ,2007 ) allowed us to deter- 
mine the star's physical radius. This result is also 
listed in Table El 

For the ^ld fit, the errors were derived via the 
reduced minimizatio n method ( Wall fc Jenkind 



20031: iPress et al.lll992[ ): the diameter fit with the 



lowest was found and the corresponding diame- 
ter was the final 0ld for the star. The errors were 
calculated by finding the diameter at + 1 on 
either side of the minimum determining 
the difference between the diameter and + 1 
diameter. In calculating the diameter errors in Ta- 
ble|31 we adjusted the estimated visibility errors to 
force the reduced to unity because when this 
is omitted, the reduced x^ is well under 1.0, in- 
dicating we are overestimating the errors in our 
calibrated visibilities. Figure [2l shows the LD di- 
ameter fit for HD 148293. Though several points 
lie outside the 1-cr line of the diameter fit, the er- 
rors in the individual visibility points overlap the 



20% change in the fj,x leads to a change in the measured 
LD diameter is less than 1% 
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Table 2 

HD 148293's Calibrated Visibilities. 



Calib 




B 




e 












HD 


MJD 


(m 


) 


(deg) 




V 


aV 




5o4Ur.iDO 


21U, 


AO 


251, 


,2 


U, 


CIO 


0, 


.077 






oil 
211, 


, (D 


253, 


,7 


U, 


A A '\ 

.441 


0, 


.068 




55407.182 


212, 


A A 

,40 


255, 


,9 


0, 


.465 


0, 


.063 




CCA AT 1 A1 

55407.191 


212, 


A C 

,95 


258, 


,2 


0, 


.439 


0, 


.061 




55407.200 


213, 


,41 


260, 


,4 


0, 


.542 


0, 


.062 




55407.208 


213, 


,77 


262, 


,7 


0, 


.504 


0, 


.057 




55407.217 


214, 


,03 


264. 


,9 


0, 


.426 


0, 


.055 


145454 


55408.275 


156, 


,09 


243. 


,2 


0, 


.759 


0, 


.072 




55408.283 


156, 


,15 


246. 


,0 


0, 


.768 


0, 


.061 




55408.291 


156, 


,18 


248, 


,7 


0, 


.752 


0, 


.035 




55408.304 


156, 


,22 


252, 


,8 


0, 


.750 


0, 


.057 




55408.312 


156, 


,24 


255, 


,7 


0, 


.765 


0, 


.080 


147321 


55408.275 


156, 


,09 


243, 


,2 


0, 


.685 


0, 


.045 




55408.283 


156, 


,15 


246, 


,0 


0, 


.728 


0, 


.052 




55408.291 


156, 


,18 


248. 


,7 


0, 


.730 


0, 


.033 




55408.304 


156, 


,22 


252. 


,8 


0, 


.763 


0, 


.052 




55408.312 


156, 


,24 


255. 


,7 


0, 


.747 


0, 


.072 



Note. — The projected baseline position angle (9) is 
calculated to be east of north. 



Table 3 

HD 148293 Stellar Parameters. 



Parameter 


Value 


Reference 


[Fe/H] 


+0.08 


Cavrel de Strobel et al. C 1997) 


V magnitude 


5.25 


Mermilliod C1991) 


K magnitude 


2.83 ± 0.11 


NeuEebauer & Leighton (1969) 


Av 


0.04 


Famaev et al. f2005) 


BC 


0.36 ± 0.10 


Alonso et al. ( 1999) 


i^BOL (10-** erg s-i cm-2) 


28.9 ± 2.8 


Calculated here 


6'uD (mas) 


1.439 ± 0.059 (4%) 


Measured here 


6'ld (mas) 


1.480 ± 0.060 (4%) 


Measured here 


J^linear {Rq) 


14.3 ± 0.6 (4%) 


Measured here 


TeS (K) 


4640 ± 100 


Brown et al. ( 1989) 


Teff (K) 


4460 ± 141 (3%) 


Measured here 
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diameter fit itself. 




BAbELINE/WAVELENC 



Fig. 2.— HD 148293 LD disk diameter fit. The 
solid line represents the theoretical visibility curve 
for a star with the best fit ^ld , the dotted lines are 
the ItT error limits of the diameter fit, the filled cir- 
cles are the calibrated visibilities, and the vertical 
lines are the measured errors. 

Once ^LD was determined interferometrically, 
the Tcff was calculated using the relation 



1 



Cff! 



(2) 



where i^soL is the bolometric flux and tr is the 
Stefan-Bolzmann constant. i^eOL was determined 
in the following way: HD 148293 's V and K 
magnitudes were dere ddened using t he ex tinc- 
tion curve described in Cardelli et al. I (Il989l) and 
its interstellar absor ption (^v) value was from 
Famaev et al. I (l2005h . The intrinsic broad-band 
color {V — K) was calculated and the bolometric 
correction (BC) was determined by interpolating 
between the [ Fe/HJ = 0.0 and -1-0.2 tables from 
Alonso et al.l ( 19991) . They point out that in the 
range of 6000 K > Tes > 4000 K, their BC calibra- 
tion is symmetrically distributed around a ±0.10 
magnitude band when compared to other calibra- 
tions, so we assigned the BC an error of 0.10. The 
bolometric flux was then determined by applying 
the BC and the Teg was calculated. See Table |3] 
for a summary of these parameters. 

4. Results and discussion 

As a check to our measured diameter, limb- 
darkened angular diameters were estimated using 



two additional methods: (1) by producing an SED 
fit (see Figure [3]) as described in Sectio n 2, w here 
UBV photometry is from M ermilliod (|l99ll ). RI 
photometry is from lMonet et al.l (l2003i): and JHK 



using the relationship described in iKervella et al 



photometry is from ICutri et al.l (12003): and (2) 



([2OO4I ) between the [V - K) color and log 6'ld. 
Our measured 6'ld is 1.480 ± 0.060 mas, the SED 
fit estimates 1.418 ± 0.085 mas, and the color- 
diameter relationship produces 1.465 ± 0.692 mas. 
Because 2MASS measurements saturate at mag- 
nitudes brighter than ~3.5 in the if-band even 
when using the shortest exposure timeQ, we used 
the K magnitude from the Two-M icron Sky Sur- 
vey ( Neugebauer fc Leighton 1969[) for the color- 
diameter determination. 




10000 20000 

WAVELEKGTH (Angstroms) 



Fig. 3.— HD 148293 SED fit. The diamonds 
are fluxes derived from UBVRIJHK photome- 
try (left to right) and the solid line is the Kurucz 
stellar model o f a star with Teg = 4571 K a.nd lo g 
g = 2.41 from lAllende Prieto fc Lambert! (|l999f ). 
There were no quoted errors for the U BVRI mea- 
surements, Gj ~ 7%, Oh ^ 7, and crx ^ 12%, 
which are not indicated on the plot. 

The main sources of errors for the three meth- 
ods are uncertainties in visibility measurements for 
the interferometric measurement, uncertainties in 
the comparison between observed fluxes and the 
model fluxes for a given Toff and log g for the 
SED estimate, and uncertainties in the K mag- 
nitude for the color-diameter determination. All 



■^Explanatory Supplement to the 2MASS All Sky 
Data Release and Extended Mission Products, 
|http : / / www, ipac .caltech.edu/2mass/releases/allsky/doc/[ 
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three diameters agree within their errors but our 
interferometric measurements provide an error ap- 
proximately 40% and 12 x smaller than the latter 
two methods, respectively. 

We also wanted to compare our Toff with those 
obtained from the literat ure (see Table H|. Val- 
ues range fr om 4390 K ( Lafrasse et al. l2010[ ) to 



4693±24 K (|Wu et al.M201l[ ). The methods are 



varied and are based on spectroscopic observa- 
tions, photometric measurements, spectral typing, 
color-Toff relationships, the correlation between 
absorption line features and stellar parameters, 
and the star's position on the color-magnitude di- 
agram. Most temperatures do not have errors in- 
dicated in the source papers. Our measurement 
of 4460±141 K is on the cooler side of the values 
from the literature, though the entire range spans 
only -300 K. 

With our newly calculated Toff, we are able 
to place HD 148293 on the H-R diagram from 
CBOO. Figure |4] shows a reproduction of their Fig- 
ure 1 with Li-rich giant stars plotted along with 
evolutionary tracks for a range of s tellar masses. 



We used the YZVAR stellar model (Bertelh et al 



[2008, 2 00i) for HD 148293 's m etallicity of +0.08 
(fCavrel de Strobel et ai]|l997t) . 

CBOO hypothesized that HD 148293 is at the 
red- giant bump, even t hough the Toff they used 
from lBrown et al.l (|l989h places the star slightly to 
the left of the bump. They claim the temperature 
shift required to place HD 148293 in the red-giant 
bump is not unreasonable (~200 K), given the un- 
certainties in the Toff measurement. Our Toff is 
180 K cooler than the value they used and places 
HD 148293 closer to and within the error bar of 
the red-giant bump, supporting their hypothesis. 
CBOO do not seem to believe that the slight shift 
in Teff significantly affects the Li abundance cal- 
culated by Brown et al., so it would appear likely 
that we are seeing HD 148293 during a very brief 
stage in its evolutionary process where Li is being 
produced or was produced very recently. 

The CHARA Array is funded by the National 
Science Foundation through NSF grant AST- 
0606958 and by Georgia State University through 
the College of Arts and Sciences, and by the W.M. 
Keck Foundation. STR acknowledges partial sup- 
port by NASA grant NNH09AK731. This research 
has made use of the SIMBAD database, operated 




Fig. 4. — H-R diagram for the Li-rich giants. The 
small circles are the stars from CBOO, the large cir- 
cle is their original placement of HD 148293, and 
the large triangle is our new placement. The evo- 
lutionary tracks are marked by the stars' masses, 
and the shaded region indicates the red giant 
bump region. The luminosity and temperature 
error bars for the Charbonnel & Balachandran 
points are from uncertainties in the Hipparcos par- 
allaxes and from the literature, respectively. Our 
error bars are based on interferometric observa- 
tions, uncertainties in Hipparcos parallaxes, and 
bolometric corrections. 

at CDS, Strasbourg, France. This publication 
makes use of data products from the Two Micron 
All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center /California Institute 
of Technology, funded by the National Aeronau- 
tics and Space Administration and the National 
Science Foundation. 
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